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INTRODUCTION 

The  objective  of  this  program  was  to  investigate  the  compressibility  of  frozen  soils  and 
ice  in  the  pressure  region  of  0  to  30  kbars  at  -10C.  The  data  were  obtained  for  use  in  model  cal¬ 
culations  for  mound  and  cavity  growth  and  shock  wave  transmission  during  a  nuclear  or  high-ex¬ 
plosive  craterinfTevent  in  frozen  ground.  The  wide  variation  of  the  properties  of  in  situ  frozen 
soils  with  location  precludes  the  direct  use  of  the  reported  test  results  for  a  particular  site. 
Therefore^  emphasis  has  been  placed  on  isolating  the  physical  properties  of  frozen  soils  that 
determine  the  compressibility  so  that  the  compressibility  of  in  situ  frozen  soils  can  be  predicted 
from  convenient  measurable  parameters. 

Frozen  soils  consist  of  a  matrix  of  mineral  grains  with  pore  spaces  that  may  be  filled  with 
ice  (water)  and  air.  The  problem  of  prediction  of  compressibilities  falls  into  two  categories: 
those  for  saturated  frozen  ground  and  those  for  partially  saturated  frozen  ground.  For  saturated 
frozen  ground  the  compressibility  can  be  predicted  from  the  constituent  components,  ice  and  min¬ 
eral.  The  required  parameters  are  the  volumetric  proportions  and  the  compressibilities  of  the  min¬ 
eral  components  and  the  ice.  For  partially  saturated  frozen  ground,  in  addition  to  these  parameters, 
the  closure  of  the  pore  spaces  due  to  the  crushing  of  the  grains  must  be  evaluated.  This  problem 
remains  partially  unsolved.  As  will  be  demonstrated,  this  appears  to  be  significant  only  for  soil 
with  a  low  degree  of  saturation. 

Bridgman  (1911,  1912,  1914,  and  1937)  was  the  first  to  investigate  the  compressibility  of 
water  and  its  many  ice  phases  to  high  pressures.  His  work  is  the  basis  for  this  discussion.  Other 
investigators  (Whalley  et  a i.,  1966;  Brown  and  Whalley,  1966,  Wilson,  et  a 1.,  1965  )  have  studied 
the  phase  diagram  of  water  in  the  pressure-temperature  plane  but  have  not  made  volume  change 
observations. 

Numerous  investigators  (Adams  and  Williamson,  1923;  Bassett  et  a I.,  1968;  Bridgman,  1964, 
Brown  et  a J.,  1967;  Stephens,  1964;  Stephens  and  Lilley,  1966;  Walsh,  1965,  a.b)  have  studied  the 
compressibility  of  minerals,  soils,  and  rocks  to  high  pressures. 

A  piston-die  device  with  which  a  uniaxial  load  is  imposed  on  a  lead-encapsulated  specimen 
was  used  for  these  tests.  The  die  restrains  the  encapsulated  specimen  from  lateral  expansion. 
Load  and  deformation  measurements  are  continuously  recorded.  The  compressibility  of  gold  is 
used  as  a  standard  to  provide  a  continuous  calibration  of  the  loading  apparatus. 

The  techniques  and  operational  procedures  were  developed  with  the  aid  of  personnel  of  the 
Lawrence  Radiation  Laboratory,  Livermore,  California. 
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Figure  1.  Encapsulated  sample  and  high  pressure  die . 


Figure  2.  Test  apparatus  in  loading  machine. 
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APPARATUS 

The  experimental  procedure  and  isothermal  compression  apparatus  are  similar  to  those 
described  by  Stephens  (1964).  The  test  samples  were  nominally  1.168  mm  in  diameter  and  2.540 
mm  long.  Each  was  encapsulated  in  lead  and  compressed  between  two  1.27-cm-diam  carbide 
pistons  while  confined  in  the  bore  of  a  die  plate  (see  Fig  1  and  2).  Hardened  steel  rings  and  plugs 
prevented  extrusion  of  the  lead. 

The  piston  was  loaded  by  a  Tinius-Olsen  Universal,  scrtw-type  testing  machine  (see  Fig.. 

3)  at  a  rate  slow  enough  to  obviate  any  compressive  heating  effects  (=  40  min/cycle).  The  tests 
were  conducted  in  a  coldroom  maintained  at  approximately  -10C.  To  obtain  a  more  precise  temp¬ 
erature  control,  a  coolant  was  circulated  through  colls  attached  to  the  die  plate.  The  coolant  was 
maintained  at  a  temperature  of  -10  t  0.5C  by  a  constant  temperature  bath.  The  die  temperature 
was  sensed  by  a  glass  bead  thermistor  and  recorded  at  intervals.  The  piston  displacement  was 
measured  by  a  linear  film  potentiometer  and  the  load  by  a  load  cell.  The  load  and  displacement 
were  recorded  on  an  X-Y  plotter.  Typical  records  are  shown  in  Figures  4  and  5. 


Figure  3.  Isothermal  compressibility  test  apparatus. 


DISPLACEMENT 
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Figure  4.  Example  of  X-Y  plots  for  saturated  frozen  soil  (2nd  compression  curve). 


Figure  5.  Example  of  X-Y  plots  for  gold  (3rd  compression  curve). 
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SAMPLE  DESCRIPTION 

The  relevant  index  properties  and  gradations  of  the  materials  tested  are  shown  in  Table  1 
and  Figure  6.  The  test  specimens  included  two  soils  West  Lebanon  glacial  till  (WLGT)  and 
Ottawa  banding  sand  (OWS).  They  were  prepared  at  various  degrees  of  saturation  (see  Table  ID. 
In  addition,  tests  were  conducted  on  clear  polycrystalline  ire. 

Ottawa  haadlag  saed 

This  material  was  well  rounded,  well  sorted  quartz  sand  with  a  median  diameter  of  100  p 
and  100%  of  the  material  falling  between  74  and  149  p.  The  material  was  selected  for  its  high 
porosity  and  granular  nature,  to  allow  sands  and  gravels  to  be  modeled.  The  nature  of  the  test 
program  restricted  the  maximum  particle  size  to  149  p.  Thj  monomineralic  structure  of  the  sand 
permitted  more  meaningful  comparison  with  previous  results  for  pure  quartz. 

Weal  Lebaaoa  glacial  till 

This  material  was  an  extremely  well  paded  till  with  particles  having  a  maximum  du'meter 
of  149  g  and  a  mean  diameter  of  36  p.  This  padation  contained  a  high  percentage  of  fine-grained 
materials  in  contrast  to  the  gradation  of  Ottawa  banding  sand.  Thin  material  was  cut  from  a 
boundary  till  and  in  its  test  form  was  classified  silt  (ML)  according  to  the  Soil  Classification 
System. 

Ice 


This  material  was  a  clear  polycrystalline  columnar  ice  having  a  maximum  grain  diameter  of 
0.6  cm,  with  the  c-axls  in  the  direction  of  compression. 


us  sto  sieve  size 
100  zoo 


Figure  6.  Gradations  of  the  toil »  tested. 
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SPECIMEN  PREPARATION 

All  soil  specimens  were  prepared  in  1.168-cm-ID  lend  capsules  approximately  3.2  cm  long, 
with  a  wall  thickness  of  0.051  cm.  The  specimen  length  was  approximately  2.4  cm.  This  provided 
approximately  0.8  cm  excess  lead  for  sealing  the  capsule. 

The  West  Lebanon  glacial  till  specimens  were  compacted  wet  by  tamping  with  a  1.15-cm- 
diam  rod.  The  Ottawa  banding  sand  specimens  were  compacted  dry  by  vibration.  The  water 
content  of  each  specimen  was  adjusted  to  the  desired  value  after  compaction.  Saturated  soil 
specimens  were  wetted  under  vacuum.  All  soil  specimens  were  frozen  rapidly  in  a  -10C  coldroom 
and  allowed  to  temper  for  at  least  24  hours. 

The  ice  specimens  were  frozen  from  de-aired  distilled  watet  in  a7.5-cm-diam  plastic  tube, 
machined  to  size,  and  placed  in  the  lead  capsules.  The  ice  specimens  were  also  tempered  for  24 
hours  at  -IOC. 


TESTING  PROCEDURE  AND  DATA  EVALUATION 

The  primary  data  were  generated  by  loading  and  unloading  the  test  samples  to  and  from  a 
pressure  of  approximately  30  kbars.  The  load  and  the  corresponding  displacement  signals  were 
recorded  continuously  (Pig.  4,  5).  the  recordings  resulted  in  a  loop.  This  hysteresis  was  primarily  the 
result  of  the  internal  friction  in  the  test  specimen,  the  friction  between  the  lead  and  the  wall 
of  the  die  bore,  uid  the  irreversible  closure  of  soil  voids.  The  techniques  used  to  correct 
for  friction  required  a  closed  loop.  Thus,  each  specimen  was  subjected  to  several  compression 
cycles  until  a  closed  loop  was  obtained. 

The  friction  was  assumed  to  be  equal  upon  loading  and  unloading;  the  average  of  the  two 
traces  represented  the  compressibility  of  the  material.  The  friction  increased  with  load.  The 
friction-load  relationship  was  assumed  to  be  unique  for  each  material  and  not  influenced  by  heat 
of  compression  or  by  structural  changes  from  one  cycle  to  the  next.  For  the  partially  saturated 
soil  samples,  the  compressibility  was  obtained  from  the  first  compression  curve  inasmuch  as  the 
subsequent  compressions  resulted  in  the  irreversible  closure  of  the  soil  voids.  The  friction 
correction  for  the  first  compression  cycle  was  obtained  from  the  closed-loop  compression  curve. 

For  the  fully  saturated  specimens  and  polycrystalline  ice,  the  friction-corrected  closed-loop  com¬ 
pression  curve  was  used  to  compute  the  sample  compressibility.  It  was  assumed  that  there  was 
no  material  loss  during  the  loading  and  unloading  cycles.  In  many  cases,  the  lead  capsule  con¬ 
taining  the  soil  specimen  ruptured,  resulting  in  a  moisture  loss.  These  tests  were  disregarded. 

After  friction  corrections  were  applied,  the  loading-unloading  curves  for  the  fully  saturated 
specimens  and  for  the  ice  still  exhibited  some  hysteresis  (see  Fig.4).  This  hysteresis  is  probably 
related  to  the  rate  effects  associated  with  the  phase  transition. 

To  correct  for  mechanical  effects  external  to  the  sample  a  differential  technique  was  used 
the  compressibility  of  the  test  specimen  was  compared  with  that  of  gold  (Bridgman,  1940; 

Stephens  and  Lilley  1967).  Every  second  sample  tested  was  gold. 

To  adjust  the  displacements  of  the  gold  samples  and  the  displacements  of  the  test  sample 
to  a  common  datum  further  correction  was  applied.  Each  test  sample  was  preloaded  to  50  lb.  This 
load  was  selected  because  it  was  large  enough  to  cause  seating  of  the  seals  but  not  so  large  as 
to  cause  significant  deformation  of  the  test  sample.  The  load  was  released,  the  piston  was  removed 
and  the  depth  to  the  upper  seal  O  (See  Fig.  7)  was  meastxed  with  a  depth  micrometer.  The  piston 
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Figure  7.  Set  correction  lor  gold. 


o  Till  Specimen  6  Gold  Specimen 

Figure  8.  Specimen  configuration  before  and  after  loading. 

was  replaced  and  the  test  continued  Dial  gauge  readings  were  observed  at  a  60-lb  load  before  and 
after  loading  For  the  gold  samples,  the  difference  between  these  readings  represented  the  residual 
displacement  required  to  fill  all  voids  between  the  seals  with  lead  (the  lead  would  flow  under 
pressure).  This  piston  displacement  is  called  the  set  5.  It  is  plotted  versus  D  for  gold  in  Figure 
7,  resulting  in  the  linear  expression 


S  0.6%  0.96SD  . 


(1) 
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The  normalizing  of  the  gold  runs  to  the  soil  or  ice  runs  is  accomplished  in  the  following 
manner.  The  set  S  corresponding  to  the  0  for  the  test  sample  is  found  from  eq  1.  The  set 
correction  Sc  required  to  adjust  the  apparent  axial  deformation  of  the  gold  Ag+m  to  the  same  datum 
as  the  test  specimen  is  equal  to  the  differenca  between  the  calculated  set  for  the  test  sample  and 
the  measured  set  for  the  gold  (see  Pig  8).  The  adjusted  axial  deformation  of  the  gold  \^m  at 
any  load  F  is  thus 


The  use  of  this  relationship  reouires  that  the  same  total  volume  of  material  be  contained  between 
the  seals  for  all  runs. 

Another  correction  was  made  to  account  for  the  expansion  of  the  die  bore  under  pressure. 
Stephens  and  LiUey  (1967)  gave  a  table  relating  the  true  pressure  P  to  the  apparent  pressure  Pa. 
Pitting  a  straight  line  to  their  data  resulted  in  the  following  relationship 


P  0.988  Pa 


(3) 


The  true  change  in  the  volume  of  the  test  specimen  Al'  is  related  to  the  apparent  change 
in  the  volume  of  the  test  specimen  AP,^.  the  apparent  change  in  the  volume  of  the  gold  specimen 
APK(a.  and  the  true  change  in  the  volume  of  the  gold  specimen  Avg  in  the  following  expression 


W. 


*  AVV 


(4) 


The  apparent  change  in  the  volume  of  the  gold  is  equal  to  the  average  value  for  the  gold  runs 
before  and  after  loading. 

Figure  8  illustrates  the  specimen  configurations  before  and  after  loading.  The  assumed 
deformation  under  load  is  cylindrical.  Th<s  actual  deformed  shape  is  probably  more  barrel-like, 
but  the  difficulties  in  the  evaluation  of  suah  a  shape  lead  us  to  the  cylindrical  approximation. 
The  apparent  change  in  the  volume  of  the  test  specimen  APalB  is  found  as  follows 


W 


1  ”Dl 

—  -  (/-  A 

I 


»♦« 


l 


(5) 


Prom  eq  3  the  die  bore  diameter  at  any  pressure  D^  can  be  related  to  the  die  bore  diameter  at 
atmospheric  pressure  by  the  following  expression 


D 


2 

c 


0.988 


(6) 


Substituting  eq  6  in  eq  5  gives; 


\V 


4  \  0.9C 


988  0.988 


/ 


(7) 
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Similarly 


AV  ^°lt  t  .  V"  \ 

**■  4  \  0.988  0.988  /  (8) 

From  Bridgman's  (1940)  work,  we  find  that  Al^  at  -10C  can  be  expressed  in  terms  of  the 
true  pressure  P  and  the  initial  volume  V®  as  follows 

\Vt  (5.650  x  10  4  P-  7.235  x  10“  7P2).  <®> 

By  substituting  1.27  cm  for  eq  4.  7,  8  and  9  reduce  to  the  following  expression  for  the  true 
change  in  volume  of  the  test  specimen  AF( 

AFg  -  1.282  |A#4t|  -  Ag+aJ  4  Fj^5.650  x  10"<P  -  7.235  x  10’ 7P2)  (10) 

The  results  of  the  tests  reported  given  in  the  relative  volume  Vj|  V®  versus  true 
pressure  P  plane.  The  relative  volume  is  found  as  follows 

F»VV°  1  \V/V°.  (11) 

SB  B  » 


where  F®  is  the  volume  of  the  test  specimen  at  atmospheric  pressure  and  -IOC 

THE  COMPRESSIBILITY  OP  ICE  AT  -IOC 

Fssential  to  tl.e  study  of  the  compressibility  of  ice  is  a  discussion  of  the  phase  diagram 
of  water.  Bridgman  (1911.  1912,  1914.  1937)  published  several  articles  which  are  still  the  main 
source  on  the  phase  diagram  of  water  and  ice.  Although  the  transition  pressures  in  the  phase 
diagram  of  water  have  been  studied  by  others  since  1937,  these  studies  have  been  concerned  with 
properties  of  water  other  than  compressibility  and  specific  censity,  the  properties  of  main  interest 
here. 

Figure  9  illustrates  the  phase  diagram  of  water  as  it  is  known  today  This  figure  shows 
that  the  following  transition  will  occur  during  the  isothermal  compressibility  of  ice  at  -10C.  At 
a  pressure  of  1  11  kbars  ice  I  will  melt  to  the  liquid  phase  (water).  If  the  pressure  is  raised 
further,  the  liquid  phase  will  freexe  to  ice  V  at  a  pressure  of  4.42  kbars.  Upon  continued  pressure 
increase,  ice  V  will  undergo  a  polymorphic  transition  to  ice  VI  at  a  pressure  of  6.25  kbars. 

Finally  ice  VI  will  transform  into  ice  VUl  at  a  pressure  of  20.8  kbars.  The  densities  and  the 
specific  volumes  of  the  various  water  phases  at  different  temperatures  and  pressures  may  be 
extracted  from  Bridgman's  data  (1911.  Tables  XXV  and  XXX,  and  1937,  Table  I). 

The  specific  volume  of  ice  1  at  0C  and  atmospheric  pressure  is  given  in  Table  III  us 
1.0900  em'/g.  It  is  assumed  that  the  coefficient  of  thermal  expansion  of  icc  I  is  small,  there¬ 
fore,  the  specific  volume  of  ice  1  at  -10C  is  also  1.0900  em'/g.  From  the  same  table  the  specific 
volume  of  ice  1  and  water  at  the  transition  point  of  1.11  kbars  and  -  IOC  can  be  interpolated  as 
1.0664  and  0.9544  em'/g  respectively.  Other  points  on  the  specific  volume-pressure  curve  are 
determined  below. 

Table  IV  gives  the  specific  vi  lume  of  water  and  ice  V  at  the  transition  points  for  various 
temperatures.  The  water-ice  1/  transition  at  -10C  occurs  at  4.42  kbars.  At  a  temperature  of  -10C 
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Figure  9.  Phase  diagram  nt  water  in  the  tempera 'ure-pressure  plane. 


Table  III.  Specific  volume  of  water  tad  Ice  I  oe  the  equilibria*  carve* 


Pressure 

(kgcm,l 

Pressure 

(kbaral  t 

temp 

it) 

Specific 
wot.  of 

water 

!  cm'.' gl 

Specific 

vot. 

change 

(cm*  g) 

Specific 
vot.  of 

ice 

ice*  gi 

0 

0.0 

0 

1.0000 

0.0900 

l  0900 

SOO 

0  49 

-  4  l 

0.9777 

.0998 

1.0778 

1000 

0.90 

-  8  7 

0.9808 

.1090 

1.0004 

1800 

1.47 

-  14.0 

0.9414 

1201 

1.0615 

2000 

190 

-  20  .  :l 

0.9283 

1318 

1.0571 

•  Bridgman  1911  tsfeea  from  Tsbte  XXX 

1  l  kg  cm1  0.98  *  lO^kbai* 

Table  IV.  Specific  voImm  of  water  aad  loo  V  oe  tbe  equilibrium  carve* 

Pressure 

Treasure 

Temp 

Specific 

velvet. 

water 

Specific 

vet. 

ckaage 

Specific 
eel.  of 
ire 

(»g/csif/ 

(kbara)  t 

(‘ti 

(cm‘/g) 

(em'g) 

(cm’  gl 

3800 

3.43 

-17.0 

0.8870 

0.0788 

0.8005 

4000 

3.92 

-13.0 

.8781 

.0733 

.8048 

4800 

4.41 

- 10. 1 

8894 

0881 

.8013 

8000 

4.90 

-  7.0 

.8810 

0034 

7970 

8800 

8.39 

-  4.2 

.8843 

.0890 

.7983 

0000 

8.88 

-  1.0 

.8478 

.0849 

.7929 

0800 

8.37 

♦  0.0 

.8418 

0810 

7902 

*  Bridgman. 

♦  1  kg  cm1 

1911.  tabea  from  Table  XXX 

0.98  \  10'1  kbar* 

i 


12 


ISOTHERMAL  COMPRESSIBILITY  OF  F  OZEN  SOIL  AND  ICE 


Table  V.  Specific  volaae  ckaa«e  oa  the  Ice  V  •  ice  VI  eqnllibrlnm  carve* 


Pressure 
(kg/  cm’) 

Pressure 

(kbata) 

Temp. 

f’t) 

Specific 
vo /. 
c Usage 
(cm'/g) 

6365 

6.2* 

-20.0 

0.03809 

6370 

6.25 

-15.0 

.03828 

6374 

6  25 

-10.0 

X  3847 

6377 

6.25 

-  5.0 

.03866 

6381 

6.25 

0.0 

.03886 

*  Bridgman  1911,  taken  from  Table  XXV 
t  Ikg'cm'  0,98  x  10“  kbars 


Table  VI.  Specific  volute  chaage  u  the  Ice  VI  -  Ice  VO  eqelllbrlia  carve* 


Specific 
vo  I. 


Pressure 

(kg/cm‘) 

Pressure 

fkbarslt 

Temp. 

r't) 

chtng* 

(cm’/g) 

22.000 

21.6 

-  0  0 

0.0667 

22.250 

21.8 

20.0 

.0570 

22,350 

21.9 

40.0 

0573 

*  Bridgman.  1937,  taken  from  Table! 
t  I  kg  cm  ’  0.98  x  10*’  kbits 


and  a  pressure  of  4.42  kbars  the  specific  volume  of  water  Is  Interpolated  as  0.8688  cmVg.  At  the 
same  temperature  and  pressure  the  specific  volume  of  ice  V  Is  interpolated  as  0.8012  cm'/g. 

From  here  on  It  becomes  more  difficult  to  find  reliable  data.  Although  Bridgman  was  most 
successful  in  obtaining  the  volume  changes  occurring  at  the  phase  changes,  he  had  great  difficulty 
measuring  the  compressibilities  of  the  various  phases  of  water.  He  did,  however,  develop  specific 
volume  data  for  Ice  and  water  along  their  equilibrium  curves. 

From  the  equilibrium  data  for  the  Ice  V  -  water  transition  (Table  IV)  it  Is  possible  to 
approximate  the  specific  volume  of  Ice  V  at  -10C.  At  a  pressure  of  6.25  kbars  and  a  temperature 
of  «0.2C  the  specific  volume  of  ice  V  is  interpolated  as  0.7907  cm’/g.  It  is  assumed  that  changes 
in  pressure  have  a  much  greater  effect  than  changes  in  temperature  on  the  specific  volume  of  ice 
V.  The  volume  change  associated  with  the  temperature  difference  of  10.2C  is  neglected  and  the 
specific  volume  of  ice  V  at  6.25  kbars  and  -10C  is  assumed  to  be  0.7907  cm’/g.  For  the  phase 
transition  of  ice  V  to  ice  VI  at  -10C,  Table  V  gives  the  change  in  volume  as  0.03847  cm'/g.  The 
specific  volume  of  ice  VI  at  -  10C  and  6.25  kbars  is  thus  0.7907-0.0385  «  0.7522  cm'/g. 

The  volume  of  ice  Vlll  at  -10C  and  20.8  kbars  is  the  next  calculation  to  be  made.  Bridgman 
did  not  recognize  this  phase  of  ice  but  thought  it  to  be  ice  Vll.  Ice  Vlll  was  first  observed  by 
Brown  and  Whalley  (1966)  and  Whalley  et  a /  (1966),  using  dielectric  techniques.  They  did  not 
provide  data  for  ice  Vlll.  However,  they  did  suggest  I'tat  thp  volume  change  associated  with  the 
ice  Vll  to  ice  Vlll  transition  is  very  small  (0  t  ■=  2.78  x  10M  cm’/g  at  the  ice  VI-VU-VHI  triple 
point) 

With  this  in  mind,  we  can  approximate  the  volume  of  ice  Vlll.  From  Bridgman's  work,  we 
find  that  ice  VII  has  a  specific  volume  of  approximately  0.60  cm'/g  at  room  temperature  and 
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Table  Vll.  Compressibility  ot  lee  * 


Phase 

Pressure 

(kbits) 

Specific 

vol. 

(cm'/g) 

Relative 

vol 

Temperature  •  10C 

Ice  I 

0.0 

1  0900 

1.000 

Ice  1 

1.11 

1.0684 

0  978 

Water 

1  11 

0  9544 

0  876 

Water 

4  42 

0.8688 

0  796 

Ice  V 

4.42 

0.8012 

0.736 

Ice  V 

6  25 

0.7907 

0.726 

Ice  VI 

6  25 

0  7522 

0  890 

Ice  VI 

20  8 

0.702 

0.644 

Ice  VIII 

20  8 

0  845 

0  592 

Ice  VIII 

49. 1 

C  80 

0  55 

Temperature  ♦  5C 

Water 

0.0 

0.9999 

1  0000 

Water 

8  87 

.8370 

0.8370 

Ice  VI 

8.87 

.7488 

0.7488 

Ice  VI 

21.7 

.702 

0.702 

Ice  VIII 

21.7 

.845 

0.645 

Ice  VIII 

49  1 

.60 

0.60 

*  Calculated  from  experimental  data  reported  by  Bridgman  (1911,  1937). 

49. 1  kbars.  We  assume  that  the  temperature  effect  on  the  specific  volume  of  ice  VII  is  small  in 
comparison  with  pressure  effects.  Furthermore,  we  assume  that  the  volume  change  associated 
with  the  phase  transition  of  ice  Vll  to  ice  V1Q  is  negligible  and  that  ice  VIM  exhibits  the  same 
compressibility  as  ice  Vll.  Thus,  the  specific  volume  of  ice  Vlll  at  49.1  kbars  and  -IOC  is  0.60 
cm'/g.  The  volume  change  associated  with  a  pressure  increase  from  19.6  to  44.1  kbars  for  ice 
Vll  is  given  by  Bridgman  (1937)  as  a  mean  value  of  0.039  cm'/g.  We  assume  that  this  value  is 
representative  in  the  pressure  range  of  20.8  to  49.1  kbars  at  -10C.  The  resulting  volume  change 
for  ice  Vlll  at  -10C  from  49. 1  to  20.8  kbars  is  0.046  cm'/g.  The  specific  volume  for  ice  Vlll  at 
-10C  and  20.8  kbars  is.  thus,  calculated  to  be  0.60  *■  0.045  0.645  cm'/g. 

The  volume  change  of  ice  Vll  to  ice  VI  at  the  transition  temperature  of  -10C  and  21.4  kbars 
is  extrapolated  from  Table  VI  to  be  0.057  cm'/g.  We  assume  that  this  value  holds  for  the  volume 
change  of  ice  Vlll  to  ice  VI  at  -10C  and  20.8  kbars.  The  resulting  specific  volume  for  ice  VI  at 
-10C  and  20.8  kbars  is  0,645  •  0.057  0.702  cm’/g.  This  completes  the  calculationo  for  the 

compressibility  of  ice  at  -10C  from  0  to  49.1  kbars.  The  results  are  tabulated  in  Table  Vll  and 
illustrated  in  Figure  10.  The  compressibility  of  each  phase  is  assumed  to  be  linear. 

Similar  calculations  can  be  performed  at  other  temperatures.  In  Figure  10,  the  compres¬ 
sibility  of  water  at  *5C  is  plotted  along  with  the  compressibility  of  ice  at  -10C.  Ice  in  more 
compressible  than  water,  therefore,  saturated  frozen  ground  would  be  expected  to  be  more  com¬ 
pressible  than  saturated  unfrozen  ground. 

Figures  11  through  13  and  Tables  Al-AUI,  Appendix  A.  give  test  results  for  lhe  isotheimal 
compressibility  of  ice  The  compressibility  as  predicted  from  Bridgman's  data  is  plotted  far 
comparison.  The  variation  of  the  test  results  from  Bridgman’s  results  can  be  attributed  primarily 
to  rate  effects  Bridgman  (1911)  reported  that  the  change  in  the  volume  of  the  phase  changes  was 
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Figure  10.  Predicted  Isothermal  compressibl'ity 
ol  ice  at  -IOC  and  water  at  * SC . 


Figure  12.  Isothermal  compressibility  of 
Ice.  specimen  43. 


Figure  II.  Isothermal  compressibility  ofico, 
specimen  II. 


Figure  13.  Isothermal  compressibility  of 
Ice.  specimen  44. 


time  dependent;  e.g.,  the  time  for  completion  of  the  liquid-solid  reaction  was  about  two  h oil's  on 
the  ice  I-llquid  boundary.  The  solid-solid  reactions  were  nearly  explosive  near  the  triple  point 
but  were  slow  at  other  points.  The  time  allotted  for  the  tests  reported  here  («  40  min/cycle)  was 
probably  Inadequate  to  allow  complete  phase  transformations  at  the  transition  pressures. 

The  slightly  greater  compressibility  of  the  rer  alts  reported  might  be  explained  by  the  pres¬ 
ence  of  microscopic  air  bubbles.  An  included  air  volume  of  approximately  1%  would  account  for 
the  differences. 

Bridgman  (1911)  found  that  ice  V  nucleated  only  in  the  presence  of  glass  splinters.  Two 
of  the  three  ice  specimens  tested  showed  the  liquid-ice  V  transition.  However,  for  specimen  43 
(Pig  12)  ice  V  did  not  nucleate  and  the  liquid  phase  froze  to  ice  VI  upon  loading  beyond  the  liquid- 
ice  V  transition.  In  all  tests  ice  V  was  observed  on  the  unloading  cycle;  this  is  consistent  with 
Bridgman's  results. 
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THE  COMPRESSIBILITY  OF  FROZEN  SAND  AND  SILT  AT  -IOC 

The  results  of  the  Isothermal  compressibility  tests  on  froren  sand  and  silt  are  given  in 
Appendix  A  and  in  Figures  14-28. 

It  was  suggested  in  the  introduction  that  only  a  few  material  properties  were  needed  to 
estimate  the  compressibility  of  boxen  ground.  This  problem  has  been  discussed  by  Brace  (1965) 
and  Stephens  (1964)  for  rocks.  The  general  approach  has  been  to  average  the  compressibilities  of 
the  minerals  maUng  up  the  rock.  Two  kinds  of  averages  have  been  employed:  the  Reuss  average 


Figure  14.  Isothermal  compressibility  of 
fully  saturated  Ottawa  bunding  sand, 

specimen  1. 


Figure  15.  Isothermal  compressibility  of 
fully  saturated  Ottawa  banding  sand, 
specimen  18. 


on  09  io 


To 

Figure  16.  Isothermal  compressibility  of 
fully  saturated  Ottawa  banding  sand, 
apecimea  84. 


Figure  17.  Isothermal  compressibility  of 
fully  saturated  Ottawa  banding  sand, 
apecimea  87. 
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Figure  18.  Isothermal  compressibility 
o!  lully  saturated  West  Lebanon 
glacial  till,  specimen  32. 


Figure  19.  Isothermal  compressibility 
ol  lully  saturated  West  Lebanon 
glacial  till,  specimen  42. 


Figure  20.  Isothermal  compressibility  ol 
partially  saturated  Ottawa  banding 
sand,  specimen  34. 


Figure  21.  Isothermal  compressibility  ol  partially 
saturated  Ottawa  banding  sand, 
specimen  17. 
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Figure  22.  Isothermal  compressibility  ol  partially 
saturated  Ottawa  banding  sand, 
specimen  28. 


Figure  23.  Isothermal  compressibility  o!  partially  saturated  Ottawa 
banding  sand,  specimen  33. 
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Figure  iS,  laolheraal  compreeelbllltg  ot  dry  Ouaw a  beading  send. 

apeclmen  SI, 
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Figure  26.  Isothermal  compressibility  ol  dry  Ottawa  banding  sand, 


°0  7  0  8  0.9  10 


l- _ i _ _ — t-  l 
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V/v0 

Figure  27.  Isothermal  compressibility 
ol  partially  saturated  West  Lebanon 
glacial  till,  specimen  16. 


V/vo 

Figure  28.  Isothermal  compressibility 
ol  partially  saturated  West  Lebanon 
glacial  till,  specimen  40. 


and  the  Voigt  average.  The  Reuss  average  B,  is  given  by: 

B,  ’  VaH  4  ‘i8*  ♦  V,  Bc  ♦  •••  02) 

where  V;.  1^.  etc.  are  volunie  percentages  of  the  minerals  in  the  rock,  and  B  .  B^.  etc.  are  the 
volume  compressibilities  of  the  minerals.  The  Reuss  average  provides  an  upper  bound.  The  lower 
bound,  the  Voigt  average  6,  is  given  by: 

1/BV  -  Vt/Bt  ♦  V«b  *  K/Bc . 


(13) 
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The  R»iuss  averaging  method  assumes  uniform  distribution  of  stress  throughout  the  matrix  and  neg¬ 
lects  the  difficulty  in  fitting  the  distorted  grains  together.  The  Voigt  model  assumes  a  uniform 
strain  throughout  the  matrix  and  neglects  the  non-equilibrium  stress  conditions  that  exist.  If 
discrepancies  occurred  between  the  predicted  compressibilities  given  by  eq  12  and  13  and  the 
measured  values,  porosity  or  alterations  of  the  minerals  in  the  rock  were  listed  as  causes. 

The  prediction  of  the  compressibility  of  frozen  soil  is  analogous  to  that  of  rock.  Fully  and 
partially  saturated  frozen  ground  are  discussed  separately. 


THE  COMPRESSIBILITY  OF  FULLY  SATURATED  FROZEN  GROUND 

A  saturated  frozen  ground  can  be  considered  as  a  rock  made  up  of  two  components-  ice  and 
mineral  particles.  The  compressibility  of  saturated  frozen  ground  can  be  estimated  by  either  eq  12  or 
13  if  we  assume  that  the  ice  flow.,  plastically.  The  problem  is  complicated  by  the  phase  tt  in- 
siticus  that  ice  undergoes  when  compressed  to  30  kbars  at  -10C.  This  fact  preclude)  the  cirect  use 
of  eq  12  or  13.  A  consequence  of  eq  12  is  that  the  total  volume  change  Al'  is  the  mere  sum  of 
the  volume  changes  in  the  component  materials,  the  mineral  particles  Al^  and  the  ice  Al{  .  l.e. 

AF,  -  AVm  .  AF,.  (14) 

With  knowledge  of  AF^  and  AF  AF§  can  thus  be  computed  at  any  pressure  and  temperature.  Val¬ 
ues  for  AFm  can  be  obtained  from  data  given  by  Brace  (1965),  Stephens  and  Lilley  (1966),  and 
Stephens  (1964). 

A  general  equation  for  the  compressibility  of  frozen  ground,  with  porosity  as  the  main  pa¬ 
rameter,  can  now  be  derived.  Porosity  n  is  defined  as  the  ratio  of  the  volume  of  voids  F°  to  the 
total  volume  V ®  at  atmospheric  pressure' 

n  -  V°/V°t.  (15) 

For  a  porous  medium  saturated  with  ice,  eq  15  becomes: 

n  -  Ff/Fj  (16) 

where  F°  is  the  volume  of  ice  at  atmospheric  pressure. 

AVj  at  pressure  P  is  given  by 

AF,  -  nV°  (1  -  vf/vf)  (17) 

where  v|Vv®  is  the  value  that  can  be  obtained  from  Table  VII,  at  any  pressure.  The  volume  of 
the  mineral  solids  F^  is: 

f®  (18) 

and  AFm  is  given  by 

AFm  =.  (F®  -  F°)(aP  -bP2) 

ro  si 


(19) 
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wliere  a  and  b  are  compressibility  coefficients  and  Pis  pressure  in  kbars.  Values  fa  a  range 
from  2.68  *  10*’  fa  quartz  to  1.01  x  10*'  fa  augite  (Brace.  1965).  Values  fa  b  range  from  24  x  HP 
for  quartz  to  3.9  *  10**  for  calcite. 

The  total  volume  change  caused  by  pressure  P  is  given  by: 

W.  nV\  (  1  -  vf/v°  )  *(F°  -  F?)(aP  -  bP*). 

The  relative  volume  ratio  Fj/F^  is  usually  plated  versus  P. 

FP/Fj  .  t\Vt/V°t  .  1  _  n  (1  -  vP/v°)  -  (1  -  n)  ♦ 

Figures  14-19  show  the  compressibilities  of  saturated  frozen  Ottawa  banding  sand  and 
West  Lebanon  glacial  till  with  the  predicted  compressibilities.  The  compressibilities  fa  ice 
obtained  from  Bridgman’s  wak  are  used  fa  the  prediction.  Other  than  in  the  regions  of  the  phase 
changes,  *he  differences  are  subtle.  The  differences  noted  at  the  phase  changes  are  of  the  same 
nature  as  those  observed  for  pure  ice.  Howeva,  the  ice  Vl-ice  VIII  transition  appears  to  begin  at 
a  higher  pressure  than  expected  (22.5  kbars  vs  20.8  kbars).  This  again  may  be  the  result  of  the 
time-dependent  behavia  of  the  reaction.  But  it  may  also  be  caused  by  a  nonhomogeneous  pressure 
distribution:  i.e.,  the  time  is  not  sufficient  fa  the  pressure  on  the  mineral  component  and  that  on 
the  ice  to  equilibrate  so  the  mineral  component  takes  a  higher  pressure. 

THE  COMPRESSIBILITY  OF  PARTIALLY  SATURATED  FROZEN  GROUND 


(20) 

Thus 

(aP  ♦  bPs).  (21 J 


By  definition,  unsaturated  frozen  ground  consists  of  a  mineral  phase,  a  gas  phase,  and  an 
ice  phase.  The  degree  of  saturation  S(  is  defined  as  the  ratio  of  the  volume  of  the  voids  filled 
witli  ice  F°  and  the  total  voids  volume  F^.  Thus 


F°  S,F®  (22) 

F° can  be  expressed  in  terms  of  the  total  volume  F^  and  the  porosity  o  by 

F°  .  S,nF°.  (23) 


In  an  analogy  with  our  analysis  fa  saturated  frozen  ground,  the  changes  in  the  volume  of 
unsaturated  frozen  ground  can  be  written  as  the  sum  of  the  true  change  in  the  volume  of  the  mineral 
material  AF,.  of  the  ice  AV'  and  of  the  air  voids  AFt>  Thus 


AF  .  AFm  ♦  AV'  ♦  AFt. 

In  a  partially  saturated  soil,  the  volume  of  the  air  voids  V°k  is  given  by 

We  will  assume  fa  now  that  all  voids  close  with  slight  pressure.  Then 
\F»  *  F®  (1  -  St)  n  F^. 


(24) 


(25) 


(26) 
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In  an  analogy  with  eq  17,  the  volume  change  of  the  ice  at  any  pressure  in  an  unsaturated 
frozen  ground  is  given  by 

s,»<  “  ■?'»?> 

and  AFm  is  given  by 

Wm  ’  {V°n  -  O  <aP  *  bP  2*> 

Hence 


(27) 

(28) 


\F  =  ( 1  S,)  n  V*  *  S(  n  F®  (1  -  v{*/v®)  ♦  O'®  )  (aP  .  bP  \  (29) 

and 

V'P/V'®  i  -  1  (1-S()  n  -  St  n  (i  v‘*  v®)  ( i  -  n)  (aP  *  bP  2).  (30> 

The  isothermal  compressibility  of  unsaturated  frozen  ground  is  thus  a  function  of  the  degree  of 
saturation  with  ice  S|  the  porosity  n.  and  the  isothermal  compressibility  of  ice  and  mineral  par¬ 
ticles. 

In  Figures  20-28,  the  results  for  the  compressibility  of  partially  saturated  Ottawa  banding 
sand  and  West  Lebanon  glacial  till  are  plotted  with  the  predicted  compressibility.  The  values 
used  for  a  and  b  in  eq  30  were  2.68  *  10'*  and  24  x  10‘‘ ;  P  was  in  kbars.  Specimens  with  various 
degrees  of  saturation  were  tested. 

For  the  partially  saturated  (S,  -  58%)  West  Lebanon  glacial  till  (Fig  27  and  28)  it  appears 
that  the  air  void  closure  occurred  at  some  pressure  below  2  kbars.  The  air  void  closure  as  ob¬ 
served  is  somewhat  complicated  and  is  partially  obscured  '>y  the  phase  change  that  occurs  at  1.11 
kbars.  The  predicted  compressibility  is  plotted  for  comparison.  At  the  higher  pressures  there 
appears  to  be  some  deviation  of  the  test  results  from  the  predicted  values.  Moreover,  on  the 
loading  cycle  the  phase  changes  nucleate  at  a  higher  pressure  than  predicted.  This  occurrence 
indicates  that  complete  air  void  closure  has  not  taken  place  or,  what  is  more  likely,  that  the 
mineral  component  is  sifcjected  to  a  greater  stress  fhan  the  ice.  This  behavior  is  undoubtedly 
dependent  upon  the  rate  of  compression. 

The  compressibilities  for  partially  saturated  Ottawa  banding  sand  are  plotted  in  Figures 
20-26.  Again  the  predicted  values  are  shown  for  comparison.  It  appears  that  the  air  void  closure 
is  not  obtained  at  2  kbars.  For  75%  and  50%  saturated  specimens  closure  does  not  appear  to  occur 
until  approximately  10  kbars.  The  phase  changes  are  not  well  defined  nor  do  they  appear  to  be 
complete.  Again,  the  phase  changes  nucleate  at  higher  pressures  than  predicted.  On  the  unloading 
curves  the  phase  changes  appear  as  predicted,  although  tb?  rate  dependence  is  still  in  evidence. 
The  difference  between  the  test  results  and  the  predicted  values  for  the  compressibilities  of  par¬ 
tially  saturated  sand  is  small  (approximately  1%  of  the  initial  volume).  This  difference  is  well 
within  the  expwimental  accuracy. 

It  can  be  observed  in  the  loading  curves  that  a  significant  pressure  is  required  to  close  the 
air  voids.  The  prediction  of  the  closure  of  the  air  voids  as  a  function  of  pressure  appears  to  be 
difficult. 

The  compressibility  of  tocks  with  small  cracks  or  with  spherical  pores  was  investigated 
by  Walsh  (1965  a.b).  He  investigated  analytically  the  elastic  behavior  of  solids  with  cracks  run¬ 
ning  through  them.  Important  parameters  are  the  shape  and  direction  of  the  crack.  Small  differences 
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in  crack  length  with  direction  could  lead  to  significant  differences  in  linear  compressibility. 
According  to  Walsh,  the  hydrostatic  pressure  necessary  to  close  an  elliptic  cavity  is 

P  =  Ea  (31) 

where  E  is  Young’s  modulus  and  a  is  the  ratio  of  miner  to  major  radius  of  the  cavity.  A  spherical 
cavity  requires  a  larger  pressure  for  closure.  Equation  31  applies  to  a  dry  homogeneous  rock,  or 
to  a  dry  nonhomogeneous  rock  if  E  is  the  Young's  modulus  of  the  weakest  material.  Equation  31 
also  assumes  that  the  cavity  is  closed  by  elastic  deformation.  If  the  crack  is  closed  by  plastic 
flow  or  by  brittle  failure,  eq  31  does  not  apply. 

On  the  other  hand,  partially  saturated  frozen  ground  consists  of  ice  and  an  unconsolidated 
mineral  matrix.  The  first  reduction  in  air  void  volume  would  occur  not  by  elastic  deformation  or 
by  plastic  flowbut  by  rearrangement  of  the  mineral  particles.  This  would  be  followed  or  accompa¬ 
nied  by  crushing  of  the  individual  particles.  The  extent  of  the  rearrangement  and  the  crushing 
would  be  governed  by  the  amount  of  ice  present.  As  the  mineral  particles  are  being  rearranged  and 
crushed,  the  ice  would  flow  plastically.  At  some  pressure  the  icr  would  completely  fill  the  voids 
and  be  subjected  to  the  same  overall  pressure  as  the  particles.  The  compressibility  then  would 
be  governed  by  the  deformation  of  the  mineral  particles  and  the  plastic  defamation  of  the  ice. 

The  pressure  at  which  the  air  voids  close  is  influenced  by  the  degree  of  saturation  with  ice. 
The  air  voids  in  soils  with  a  high  degree  of  saturation  would  close  at  a  relatively  low  pressure 
while  those  with  a  low  degree  of  saturation  would  require  a  higher  pressure  fa  air  void  closure. 
Maeover,  the  rearrangement  and  crushing  may  be  expected  to  be  more  efficient  fa  a  wet  soil  than 
fa  a  dry  soil. 

The  Ottawa  sand  specimen  33  (Fig.  23),  with  a  25%  saturation,  exhibits  a  somewhat  different 
behavia  from  that  of  the  specimens  with  higher  degrees  of  saturation.  Closure  does  not  appear  to 
occur  until  a  pressure  of  approximately  20  kbars  has  been  reached.  Upon  unloading,  the  test  data 
follow  the  predicted  data  down  to  a  pressure  of  approximately  1  kbar.  The  phase  changes  are  not 
well  defined.  This  would  be  expected  because  of  the  small  volume  of  ice  present.  However,  at 
approximately  1  kbar,  the  test  data  leave  the  predicted  values  and  the  specimen  appears  to  expand 
elastically.  Figures  25  and  26  show  that  the  dry  Ottawa  sand  specimens  have  the  same  elastic 
expansion  fa  the  low-pressure  release  curve.  In  fact,  the  dry  and  25%  saturated  specimens  release 
to  approximately  the  same  relative  volume  (0.78).  This  relative  volume  is  the  smallest  volume 
regardless  of  the  degree  of  saturation  to  which  Ottawa  banding  sand  can  be  compressed  under  the 
test  conditions. 

The  Ottawa  banding  sand  compressibility  curves  look  much  like  those  reported  for  dry 
Monterey  sand  by  Stephens  and  Lilley  (1966).  Complete  closure  of  the  voids  is  not  obtained  and 
the  mineral  particles  defam  elastically. 

To  test  the  validity  of  the  methods  used  in  reducing  the  raw  data  fa  the  partially  saturated 
soils,  the  first  and  second  compressibilities  were  calculated  fa  a  fully  saturated  sand.  The  results 
of  this  comparison  (Fig.  29)  indicate  that  the  methods  used  to  reduce  the  data  from  the  first  com¬ 
pression  are  valid.  However,  the  phase  transition  of  ice  I  to  water  is  somewhat  obscurer  and  the 
subsequent  phase  transitions  ocor  at  somewhat  high  pressures  than  expected. 

The  results  of  three  compression  cycles  on  dry  Ottawa  banding  sand  are  illustrated  in 
Figure  30.  This  figure  shows  that  the  voids  undergo  maximum  closure  during  the  first  compression. 
The  first  release  curve  and  the  additional  compression  cycles  follow  nearly  the  same  path.  Thus, 
it  appears  that  the  working  of  the  mineral  particles  and  the  subsequent  breakdown  do  not  result  in 
a  further  change  in  compressibility. 
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Figure  29.  First  and  second  compressibilities 
of  fully  saturated  Ottawa  banding  sand, 
specimen  24. 


Figure  30.  First,  second  and  third 
compressibilities  of  dry  Ottawa 
banding  sand,  specimen  35. 
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The  isothermal  compressibility  of  indium  was  determined  as  an  additional  check  on  the 
accuracy  of  our  procedures  and  equipment.  The  procedures  followed  were  identical  to  those  used 
far  the  compressibility  of  the  frozen  sand,  frozen  silt,  and  ice.  The  compressibility  of  indium  was 
compared  with  that  of  gold.  The  results  are  given  in  Figure  31.  Results  obtained  by  Bridgman 
(1935)  and  Stephens  (1967)  are  shown  for  comparison. 


CONCLUSIONS 

The  compressibility  of  frozen  soil  is  readily  predicted  from  the  knowledge  of  material  prop, 
erties  such  as  the  degree  of  saturation  with  ice,  the  porosity,  and  the  compressibilities  of  the 
ice  and  mineral  components.  The  behavior  of  the  phase  transitions  for  partially  saturated  frozen 
soils  is  somewhat  obscured  by  the  rearrangement  and  crushing  of  tho  mineral  particles.  Materials 
with  no  ice  or  a  low  degree  of  saturation  demonstrate  elastic  rebound  on  the  release  leg  of  the 
compression  cycle.  Apparently  there  is  a  maximum  dry  density  to  which  a  particular  soil  can  be 
compressed,  and  this  density  is  somewhat  less  than  that  of  the  voidless  mineral  parent. 

The  compre8sitdlily  of  ice  is  as  predicted  by  Bridgman  (1911,  1937).  However,  lime-depend¬ 
ent  behavior  is  demonstrated  for  the  phase  transitions.  Complete  phase  transitions  were  not 
observed  at  constant  pressure. 

The  compressibilities  in  the  low-pressure  region  (below  2  kbars)  were  not  well  defined. 

This  is  primarily  a  result  of  the  test  method.  A  technique  employing  a  liquid  cell  in  the  low- 
pressure  region  is  being  pursued  and  the  results  will  be  roported  later. 
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Table  AI.  Isothermal  compressibility  of  poly-  Table  All.  Isothermal  compressibility  of  Table  AID.  Isothermal  compresaibility  of  poly¬ 
crystalline  ice  specimen  11.  polycrystalline  ice  specimen  43.  crystalline  ice  specimen  44. 

Pressure  Volume  Relative 
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Tibi*  A IV.  Isothermal  compressibility  of 
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The  Isothermal  compressibilities  of  Ice  and  partially  and  fully  saturated  sand  and 
silt  at  -1(5* C  are  presented.  The  tests  employ  a  plston-dle  device  vlth  which  a 
uniaxial  load  Is  Imposed  on  a  lead  encapsulated  specimen,  resulting  In  the  hydro¬ 
static  compression  of  the  test  specimen.  Pressures  to  30  kbars  are  obtained.  The 
compressibility  of  Ice  Is  as  reported  by  P.W.  Brldpnan.  The  various  phase  trans¬ 
formations  of  Ice  I  to  water  to  Ice  V  to  Ice  VI  to  Ice  VUI  appear  as  expected.  It 
la  shown  that  the  compressibility  of  frozen  soil  can  be  readily  predicted  from  the 
knowledge  of  material  properties  such  as  degree  of  saturation  with  Ice,  porosity, 
and  the  compressibilities  of  the  ice  and  mineral  components. 
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